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ABSTRACT 

We analyze the relationship between maximum cluster mass, M max , and surface densities of 
total gas (S gas ), molecular gas (Eh 2 ) and star formation rate (Esfr) in the flocculent galaxy 
M 33, using published gas data and a catalog of more than 600 young star clusters in its disk. 
By comparing the radial distributions of gas and most massive cluster masses, we find that 

M max oc Sg a 7 s ±0 4 , M max oc Sh2 ±01 ' and Mmax K S sfr 01 - We rule out that these correlations 
result from the size of sample; hence, the change of the maximum cluster mass must be due to 
physical causes. 



Subject headings: galaxies: star clusters 
galaxies: individual (M 33, NGC 598) 



galaxies: ISM — galaxies: spirals — stars: formation 



1. Introduction 

A long-standing problem of galaxy formation 
and evolution has been understanding the relation 
between gas surface density (S gas ) and star forma- 
tion rate (SFR; i.e., the star formation law). In 
the last half century, significant efforts have been 
undertaken to clarify this subject. Empirical cor- 
relations have been found between the disk surface 
density of star formation, Esfr, an d the gas sur- 
face density, E gas -either total, neutral or molecu- 
lar (usually base d on CO) . Examples o f these are 



J SFR 



mi), and E SF r 



AT,^ S (ISchmidd I1959L 119631 iKennicutt 



A'Y O 



with 57 gas the av- 



erage angular velocity of th e gas in the disk (e.g., 
ISilldll997t lKennicuttf[l998l ). Data are often com- 
patible with more than one correlation at a time, 
a situation that does not help to clarify what pro- 
cesses really drive star formation rates (e.g., tur- 
bulence, large scale shocks, gravitational instabil- 
ities, shear, pressure). 

The SFR is often measured through the emis- 
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sion in the Ha line, that is due to the reprocessed 
ionizing photons produced by O or early B-type 
stars, or via the non-ionizing far ultraviolet (FUV) 
flux dominated by B-type stars. In order to infer 
the SFR in this way, then, the relation between 
the mass in massive stars and the rest, i.e., the 
initial mass function (IMF), has to be known0 

Different workers ( e.g.. ISullivan et al.l l2000t 



Bell fc Kennicuttll200ltlMeurer et al.ll2009l : lLee etal 



20091 iBoselli et al.ll2009l ) have found that, on av- 
erage, for large samples of galaxies, the SFRs 
inferred, respectively, from the Ha line and the 
FUV flux are not consistent with a universal IMF. 
Instead, the ratio between Ha and FUV seems 
to decline with galaxy lum inosity or g alaxy mass 



Lee et al.l (|2009h 



Con- 



or SFR or SFR per area, 
elude that none of the following factors can cause 
this trend if acting alone: uncertainties in stel- 
lar evolution and atmospheres, effects of different 
metallicities, variations of star formation histo- 
ries, photon leakage, extinction, and stochasticity 
of massive star formation. 

One particular case of the mismatch between 



1 The assumption also has to be made that the star formation 
activity is constant during the lifetime of the massive stars 
(si 10 7 yr for O-stars, and ss 10 s yr, for B-stars). 
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the Ha and the FUV emissions occurs in the outer 
regions of disk galaxies. Historically, the Ha cut- 
off there led to the concept of a gas surface density 
threshold, below which star formation activity was 
inhibited, and the relationships shown above be - 
tween £ gas and SFR broke down (|Kennicuttlll998» . 
Recently available UV data, however, show that 
there is no corresponding "FUV cut-off' , and that 
star formation also goes on in outer spiral disks 
(jBoissier et al.ll2007t ). 

Pflamm-Altenburg fe Kroupal (|2008f ) propose 



the local integrated galaxial IMF [(L)IGIMF] the- 
ory, that is able to explain both the Ha cut-off 
in galaxy disks, and the trend in the ratio of Ha 
to FUV fluxes with galaxy mass. The (L)IGIMF 
theory deals with surface densities (e.g., the lo- 
cal IMF and SFR densi ties), and is an outgrowth 
of the IGIMF theory dKroupa fc Weidnerl l2003t 
Weidner fc Kroupal2005t Pflamm-Altenburg et al. 



20091 ). The IGIMF theory is based on the simple 
concept that star formation, in the form of em- 
bedded clusters, occurs in molecular cloud cores 



(e.g.. lLada fc Ladall2003l ). Although most of these 
embedded clusters do not survive as bound star 
clusters the expulsion of their residual gas, the 
determination of the stellar IMF of the whole 
galaxy or of a part of it at a given time reduces 
to adding the stellar IMFs of all newly formed 
embedded clusters. On the other hand, the em- 
bedded cluster maximum mass, M ec i,tna.y, seem s 



to be a function of the total SFR (Weidner et al 



20041 ). and the stellar upper mass limit of each 
cluster's IMF is a function of the total star cluster 





Elmegreen 1983 


: Weidner & Kroupa 


2006; 


Weidner et al. 


2010|). The result- 



ing IGIMF, then, depends on the SFR, which 
itself depends on the gas density. In order to 
accou nt for the Ha cut-off in exponent ial gas 
disks, Pflamm-Altenburg fc Kroupal ( 20081 ) adopt 



the ansatz M. 



ccLmax 



cx £ 



3/2 



Another correlation between maximum cluster 
m ass and SFR surface density has been proposed 
bv lBillett et alj (2002), based on the star forma- 
tion law, and assuming the ambient interstellar 
medium and the cluster-forming cloud cores are in 
pressure equilibrium. In this scenario, M ec i )lnax oc 
Sg FR , with 2/3 < r\ < 2. The first value is ex- 
pecte d in the case of eq ual (volume) density clus- 
ters ( Billett et al.l I2002T ) . while the second would 



is supported b y the very weak birth ra dius-cluster 
mass relation (M arks fc Kroupallioil . 



In view of the limited availability of sufficiently 
accurate cluster mass determinations a decade 
ago, the relationship between maximum cluster 
lumin osity and Ssfr was investigated by lLarsen 
(2002) using HST data of 6 spiral galaxies. He con- 
cluded, however, that the effects of random sam- 
pling statistics in determining the brightest ob- 
served cluster luminosities would make it hard to 
unveil the connection between physical processes 
and maximum star cluster mass, and that a study 
performed with masses, not luminosities, would 
be much better to establish how maximum cluster 
mass depends on galaxy properties. At present, 
the widely held understanding is that the most 
massive object in an ensemble of clusters scales 
with the sample size as a result of statistical varia- 
tions. On the other hand, if there is a physical link 
between the most massive clusters and the SFR, 
then a potentially very powerful new method of 
measuring star formation histories of ga laxies is 
opened up ( Maschberger fc Kroupal 2007 ). 

Here, we test the environment dependent 
IGIMF ansatz versus the stochastic sampling 
ansatz. We make a direct comparison between 
cluster mass and gas surface density in M 33, 
using data from the literature, with the aim of in- 
vestigating whether £ gas has a role in determining 
maximum cluster mass. We will from now on refer 
to maximum cluster mass as M max exclusively. 

2. Star cluster data. 



ISharma et al. ( 2011 ) have recently published 
a study of 648 clusters younger than 10 s yr in 
the disk of M 33, out to w 16 kpc of the galaxy 
center, selecte d from the Spitze r 24 ^m image 



of the galaxy (|Verlev et al.l 120071 ). ISharma et al 
obtain galactocentric radii for the clusters using 
the w arped disk model by ICorbelli fc Schneiderl 
( 19971 ). They derive ages, masses, and extinction 
corrections from the comparison between cluster 
spectra l energy distributions (SEDs) and Sta r- 



bust^ dLeitherer et al.ll999l 



20051 ) models; ISharma et al 



Vazquez fc Leitherer 



assume a distance 



to M 33 of 840 kpc (|Freedman et al|[l99 1)FI The 



occur if clusters have equal sizes (|Larsenll2002l ) , as 



2 At this distance, 1" =4.1 pc; R25, the galactocentric 
distance of the isophote with surface brightness in the B- 
band /i S = 25 mag 1 , is 8.6 kpc or 35'.4. 
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SEDs are built with aperture photometry based on 
the Galaxy Evolution Explorer (GALEX) sat ellite 
far- and near- UV data dGil de Paz et al.ll2007l): 



Ha m ap (|Greenawalt et al.ll998tlHoopes fc Walterbosl 
2000); a nd mid- and far-IR i mages (8, 24 , 70, and 



160 ^m; IVerlev et all 120071 ). ISharma et all (|201lh 

compare the masses and the 24 /mi semi-major 
axes r of their clusters; they find a large scatter 
for r < 10 pc (about 13% of their sample), and 
-Mciust cx r 2 09 ± 001 for clusters with r > 10 pc. 
We extracted the masses and galactocentric 



radii o f the clusters from Fig. 13 in Sharma et al 



(2011) using the tool Dexter (jDemleitner et al 
200ll ). Owing to different surface brightness and 



crowding characteristics, the completeness limit 
varies with radius, and goes from 800-1000 M 
within about 0.6 R25 to w 300 Af© beyond i?25- In 
what follows, we will restrict ourselves to working 
with the 258 clusters with at least 10 3 M Q . 

3. ISM data. 



Corbelli (2003) and lHever et all (j2004l ) publish, 



respectively, neutral and molecular gas radial pro- 
files of M 33, assuming a distance of 840 kpc and 
a location of the galaxy center at RA (J2000) = 
01:33:50.89 and Dec(J2000) = 30:39:36.7. In both 
works, the data are d e proje cted with the model of 
Corbelli fc Schneiderl (jl997l ). While the published 



HI profile extends as far as 1.5 R25, the H 2 one 
reaches only to 0.8 i?25- 

For the CO to H2 conversion, Hever et al.l take 
a consta nt factor X = 3 x 10 20 cm -2 (K km 
s- 1 )- 1 !! iHever et al.1 l|2004h also derive a radial 
profile of Esfr from the FIR luminosity, using 
IRAS HiRes 60 and 100 fan images of M 330 They 
find S SFR = (3.2 ±0.2) E^™ lar , but E SFR = 



(0.0035 ± 0.066) S 



3.3±0.07 



They also ascribe 



gas, total ' 

the steep slope of the correlation with total gas 
to the very shallow radial distr ibution of atomic 



gas. Once again, we use Dexter (jDemleitner et al 



20011 ) to extract the profiles; we calculate total gas 
mass surface density as E gas = 1.36(Ejj 2 + Em), 
in order to include helium. 



3 Sjj 2 = 2 m^N (H2)', the molecular hydrogen column den- 
sity is iV(H 2 ) = X /T mb (CO(l - 0))dv, where T mb is the 
main beam brightness temperature. 

4 E SFR = 3.8 X 10- 16 (2.58</ 60 >n + (/ioo}n)M Q pc- 2 yr-\ 
with {/eo)f2 an d (^ioo)n the mean 60 and 100 fim intensi- 
ties, respectively, within a solid angle f2 . 



4. Analysis and discussion 

Results are shown in Figure [1] In the left 
panel, we plot logio of cluster mass vs. galactocen- 
tric radius in units of R25, for the objects in the 
Sharma et al.l (|2011l ) sample with mass Af c i ust > 
We constructed 7 bins covering the range 
with both neutral and molecular gas data, i.e., be- 
tween the galactic center and 0.8 R25, each one 
4f04 (1 kpc, w 0.11 R 25 ) wide. The mean of 
the five most massive clusters in each bin, and 
a weighted^] linear fit to the mean are shown, re- 
spectively, as blue crosses and a blue short-dashed 
line; the error bars are those of the mean. The me- 
dian of the five most massive clusters in the same 
bins (or the third most massive cluster in each 
bin), and a linear fit to the median are displayed 
as red filled circles and a solid line; the error bars 
represent the interquartile range. Since the uncer- 
tainty in the measurement of any individual mass 
is typically a factor of 2-3, fits to the mean and 
the median of the 5 most massive clusters in each 
bin should in general be more robust indicators 
of any existing trend (in every bin there are al- 
ways more than 5 clusters more massive than the 
completeness limit). Using this mean or median is 
even more necessary in this case, since we intend to 
compare cluster masses, not with the gas surface 
densities at their position, but with azimuthal av- 
erages at their galactocentric distance. The RMS 
azimuthal variations are always below 5% of the 
CO and mid-IR emission, and below 2% of the HI 
and HI+CO emission. 

From the weighted fit to the mean, we obtain: 

log 10 Af mcan /M = (-2.2±0.2)i?/i? 25 +(5.2±0.1); 

(1) 

the fit to the median yields: 

logio M 3ld /M Q = (-2.2±0.2)i?/i? 25 + (5.1±0.1). 

(2) 

We note that there is an extremely massive 
cluster at m 0.85 R25, almost at the edge of the 
available H2 data and beyond the published radial 
profiles of Ejj2 and Esfr- This cluster, the most 



5 Three clusters at galactocentric distances 1.07 R25, 1-20 
-R25, and 1.28 -R25, with masses, respectively, 1.6 xlO 3 Mq, 
5 X 10 4 M Q , and 3.2 X 10 3 Af Q are out of the figure, to the 
right. 

6 The mean of each bin i is weighted by Wi = l/o" 2 , where 
Oi is its dispersion. 
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Fig. 1. — M 33, radial distributions. Left panel: log cluster mass. Grey dots: cluster data; blue crosses and 
blue short-dashed line: average of 5 most massive clusters in bins 4.'04 wide (i.e., 1 kpc, rs O.H-R25), and 
weighted fit; red filled circles and red solid line: median of 5 most massive clusters (i.e., effectively, the third 
most massive cluster) in the same bins, and fit. Right panel: log surface densities (thick lines) and linear fits 
(thin lines). Solid: H2; dotted: HI; dashed: total gas; solid triangles and dotted-long- dashed line: SFR. See 
text for the fitted relation. The x- and y- axes of the two panels have the same dynamic ranges, so that the 
slopes of the fits to the cluster and gas data are directly comparable. 
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massive in the galaxy with M c i ust > 10 6 M Q , sits 
over one of the brightest spots in the CO J=l-0 
map, while there is very little emission detected 
from the rest of the galaxy at the sam e galacto- 
centric distance (see iHever et al.l 12004 their Fig- 
ure 2.). The region is also quite bright at 21 
cm, and is the brightest location between 24 and 
170 urn dHever et alJ 120041 Iffippelein et~aT1 l2003t 
Tabatabaei et al.l l2007h ~ Hence, the azimuthally 



in bins, each one containing an equal number of 
objects from the subsample of 258 clusters with 



M c i ust > 1O 3 M detected bv lSharma et all (|201ll ) 



averaged emission of the corresponding annulus 
would not be representative of the star-forming 
conditions at the location of the cluster. This clus- 
ter has, instead, formed in a local major instability 
in the interstellar medium. 



The SFR and gas surface densities (jHever et al 



2004) , and our linear fits to the data are shown in 
the right panel of Figure [TJ The x-axes of both 
panels in the figure are the same, whereas the y- 
axes have the same dynamic ranges, so that the 
slopes of fits to gas surface densities and cluster 
masses are directly comparable. 
The fits to the available gas are: 

log 10 E gas = (-0.47 ± 0.01)i?/i? 25 + (1.363 ± 0.005) 
log 10 E H2 = (-1.63 ± 0M)R/R 25 + (0.84 ± 0.01) 
log 10 E SF R - (-2.17 ± 0.02)i?/i? 25 + (1.62 ± 0.01) 
log 10 E Hl = (-0.23 ± 0.02)^/^25 + (1.05 ± 0.01)13 

For M 33, then, both the mean and the median of 
the 5 most massive clusters yield: M mcanj 3 r d oc 

V4.7±0.4. 7i/r v 1.3±0.1. M ' 

^gas i JM mean,3ni °C Zj H2 , AZ me an,3rd OC 

v 1.0±0.1 

SFR 

We discard here the possibility that the change 
in mean and median maximum cluster mass with 
radius is a statistical effect, due to the number of 
clusters in each equal sized bin decreasing with 
increasing radius. If clusters are drawn purely 
randomly, or stochastically, from the same mass 
distribution function that declines with mass, and 
that has a constant upper truncation mass M u > 
10 6 Af Q , the probability of picking a massive clus- 
ter decreases with the size of the sample, i.e., with 
the number of clusters in the sample. We demon- 
strate in the following that the change in maxi- 
mum mass is due to physical causes instead. 

Figure [3 shows fits to the median mass clusters 



7 Heyer et al. find £ gas oc RjR^? , but they do not show the 
HI data they used. We have obtained the fits to the gas 
profiles and to the cluster masses in a consistent fashion. 



in M 33; the total number of bins increases from 
3 (upper left panel) to 6 (lower left panel). The 
number of clusters in each bin is indicated, and the 
maximum and median masses are shown, respec- 
tively, with black empty triangles and red filled 
circles. 

Two fits to the medians are performed: one in- 
cluding all the bins in each panel (red dashed line) , 
and one omitting the last bin (blue solid line); this 
bin includes a radial range for which there is no 
molecular gas data, and always contains the most 
massive cluster in M 33 (see above), an object 
whose characteristics do not seem to correlate with 
the average star-forming conditions at its galacto- 
centric radius. The fits have the form: 

logio M clust /M & = {(3' ± a 0/ )R/R 25 + (a' ± a a ,), 

(3) 

and their coefficients and their uncertainties are 
listed, respectively, in the top and bottom of Table 
1. When the last bin is included, the slopes (3' of 
these fits get slowly steeper, as the number of bins 
increases (see top of Table 1) and the relative im- 
portance of the last bin progressively diminishes. 
If, on the other hand, we exclude the last bin, we 
find that the slope of the fits to the median clus- 
ter mass is fi' — —2.0 ±0.3 (see Table 1, bottom), 
regardless of the number of bins. 

If, instead of the median, we consider the 90 th 
percentile cluster (i.e., the 9 th , 6 th , 5 th , and 4 th 
most massive cluster, respectively, for 3, 4, 5, and 
6 bins), the correlation (including the last bin) 
is consistent in all cases with logio Mqq%\l/Mq — 
1.2±0.2i?/i?25 + (4.7±0.1). Thus, the relation be- 
tween cluster mass and radius is slightly shallower 
than for .Mard, but still quite robust and definitely 
not flat. We note here, though, that for large sam- 
ples (for example, of a few hundred clusters) the 
mass of the 90 th percentile point might already be 
one or two orders of magnitude below that of the 
most massive clusters, and hence would not be a 
good estimator for this particular problem. 

In order to convince the reader that the intrin- 
sic mass distribution of clusters changes as a func- 
tion of radius, we perform a Kolmogorov-Smirnov 
(K-S) test on the 6 different radial subsamples pre- 
sented in the bottom left panel of Figure 2. Figure 
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Fig. 2. — M 33, log cluster mass vs. radius, fits to bins with equal numbers of clusters. The number of bins 
increases from 3 (upper left panel) to 6 (lower left panel) ; numbers of clusters in each bin are indicated. The 
bar at the bottom of each panel shows radial ranges of the bins. Empty triangles: most massive cluster in 
each bin; filled red circles: median of five most massive clusters (i.e., third most massive cluster) in each 
bin; red dashed line: fit to medians in all bins; blue solid line: fit to medians excluding last bin. Fits are 
performed only for three and more bins. 
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Table 1 

Fits to log M 3rd/M0 vs. R/R25, Figure 2 and eq. 3 



N b 




/3 ± api a' ± a a t 


RMS 


3 


86 


(-0.5 ±0.2) (5.01 ±0.08) 


0.04 


4 


64-65 


(-0.5 ±0.6) (4.9 ±0.2) 


0.2 


5 


51-52 


(-0.6 ±0.6) (4.8 ±0.2) 


0.2 


6 


43 


(-0.7 ±0.6) (4.8 ±0.2) 


0.2 


Without last bin 


3 


64-65 


(-2.00 ±0.02) (5.17 ±0.01) 


0.003 


4 


51-52 


(-2.0 ±0.3) (5.1 ±0.1) 


0.05 


5 


43 


(-2.0 ±0.3) (5.1 ±0.1) 


0.1 



Note. — Col. (1): number of bins. Col. (2): number of clusters in each bin. Col. (3): best-fit slope. Col. 
(4): best-fit intercept. Col. (5): best-fit RMS residual. 



3 shows the cumulative probability distributions 
and median radii of each subsample, and the D 
and P values for every bin pair are given in Table 
2 (higher bin number indicates larger radius). The 
visual impression, that the inner bins arc different 
from the outer ones, is confirmed by the K-S statis- 
tic: assuming that sample pairs with P < 0.05 are 
taken from different distribution functions with 
high significance, we conclude that bins 1 and 2, 
within R = 0.22 R25, are different from bins 5 
and 6, beyond R = 0.42 i?25- Bin 3 also differs 
from bin 6, such that the change of the cluster 
mass distribution function appears to be gradual, 
rather than stepwise. 

These results rule out the size of sample effect, 
according to which no correlation (/?' = 0) is ex- 
pected between maximum cluster mass and galac- 
tic radius. From the ratio <7p>/(3', it follows that 
the falsification of stochastic sampling is at a con- 
fidence level of at least 7a. 

5. Conclusions 

We have analyzed the relationship between 
maximum cluster mass and gas surface density 
in M 33, in order to explore the suggestion that 
maximum cluster mass is determined by physical 
processes, e.g., the equilibrium pressure between 
cluster formi ng cores and the ambient inte rstcl- 
lar medium (|LarserJl20o3 iBillett et al.ll2002l) . and 
since the existence of such a relationship can rec- 



oncile, via the IGIMF theory, SFR measurements 
derived, respectively, from Ha and FUV emission 
in galaxy disks. 

To thi s end, we hav e used published g as data 
of M 33 ( Corbellil[2003l iHever et al.ll2004h . and a 
catalog of more than 600 you ng star clusters in it s 
disk, also from the literature ( Sharma et al.|[2011 ) . 
Because often the most massive cluster in a range 
of galactocentric distances is formed under condi- 
tions that are not average for the annulus, we find 
that it is best, in the present annular averaging 
approach, to use the median of the five most mas- 
sive clusters (i.e, the third most massive cluster) 
in each bin for this kind of analysis. 

We have compared radial distributions, and 



-a.3±o.i 



while M; 



3rd OC 



have found that 

^SFR° 1 (i R a ran S e consistent with the expec- 
tations from pressure equilibrium considerations) . 
On the other hand, M max cx T,^ S ±0A , steeper than 
needed to explain the Ha cut-off in galaxy disks. 
Both this correlation and the steeper than average 
star formation law might be related to the shallow- 
ness of the HI profileo 



3 It is very likely that M 33 has interacted with M 31 in the 
past. The radial HI density profile of M 33 may have been 
changed away from an exponential disk as a result, and 
hence may not reflect at present the physically relevant 
conditions for star formation in a virialized self-regulated 
galactic disk in equilibrium. Another example of a transi- 
tory state of the HI gas is the Magellanic Stream, parts of 
which will most likely be re-accreted onto the Large Mag- 
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Fig. 3. — K-S test; cumulative probability distributions of mass for clusters with M c i ust > 10 3 Mq in the 
six annuli whose median radii are indicated. Red solid line: bin 1, R = O.OGIq'q^ R25', green dotted line: 
bin 2, R = O.2I+Q Q5 R25', blue short dashed line: bin 3, R = 0.261q'q3 R25', cyan long dashed line: bin 4, 
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Table 2 
k-s test d and p values 





D P 


D P 


D P 


D P 


D P 


Bin 


2 


3 


4 


5 


6 


1 


0.56 0.9168 


0.78 0.5799 


1.22 0.1005 


1.45 0.0307 


1.63 0.0096 


2 




0.67 0.7651 


1.33 0.0569 


1.45 0.0307 


1.72 0.0052 


3 






1.00 0.2694 


1.11 0.1687 


1.63 0.0096 


4 








0.56 0.9168 


1.10 0.1797 


5 










0.69 0.7271 







Note. — D and P values for bin pairs. The cell in the intersection of a row and a column contains the D and P 
parameters, respectively, of the comparison between the two bins indicated in the corresponding row and column. If 
P < 0.05, the null hypothesis that the clusters in the two bins are taken from the same mass distribution function is 
rejected. 



In order to test whether the trend of M max 
(or its proxy, M^), with galactocentric radius is 
consistent with random sampling from the cluster 
mass function, we have also measured the radial 
distribution of maximum mass in 3 to 6 bins with 
an equal number of clusters in each bin. After ac- 
counting for the presence of the most massive clus- 
ter in the galaxy at ~ 0.85i?25, whose formation 
environment is certainly not represented by the 
average conditions across the galaxy at its galac- 
tocentric distance, and for the lack of gas data 
beyond this same radius, we find exactly the same 
results as before, regardless of the width of the 
bins. A K-S test on the mass distributions in 
these bins suggests that the two bins closest to 
the galaxy center are different from the two most 
external ones, i.e., that the mass distribution func- 
tion changes with radius. 

The significant decrease of log M max with ra- 
dial distance in M 33, as a power law with in- 
dex /3' w —2.0 ± 0.3, and despite there being 
the same number of clusters per radial bin, rules 
out random sampling with extremely high con- 
fidence. This one galaxy, thus, falsifies the hy- 
pothesis that the most massive cluster masses 
scale with the size of the sample. Instead, the 
range of star cluster masses is driven by envi- 
ronmental physics. Indeed, the available data 
for M 33 suggest logi M 3rd oc log 10 S H2 oc 



ellanic Cloud once it orbits to a larger distance from the 
MW. 



log 10 £sFR- This, however, may be merely a 
(trivial) self-consistent result because, after all, 
stars may form from molecular clouds in a free- 
fall timescale and subs equently destroy the clouds 
([Hartmann et al.ll200ll ). The non-trivial challenge 
remains to answer how the interstellar medium, 
and thus mostly the HI gas, arranges itself to form 
molecular clouds. 

We thank Mark Heyer, who provided us with 
measurements of the azimuthal variations in the 
gas data, and the anonymous referee, for his/her 
positive and helpful feedback. RAGL acknowl- 
edges support from DGAPA, UNAM. 
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